Self Gravitating Instability of a Compressible Fluid Cylinder 

Ahmed E. Radwan and Gamal G. Nashed 
Mathematics Department, Faculty of Science, Ain Shams University, Cairo, Egypt 



Abstract 

The self-gravitating instability of a compressible-inviscid fluid cylinder im- 
mersed into a self-gravitating tenuous medium of negligible motion is developed. 
The stability criterion is derived based on the linear perturbation technique. Some 
previous reported works are recovered. The effect of different factors on the fluid 
cylinder instability is discussed. The compressibility has a tendency for a sta- 
bilizing the model in particular as the sound speed a is very large in the stable 
domains 1.0233928 < s < oo but comparatively small in the unstable domains 
< X < 1.0233928 where x = kR^ is the dimensionless longitudinal wavenumber 
with k is the axial wavenumber and is the radius of the cylinder. In the absence 
of the compressibility factor the unstable domain is found to be < x < 1.0678. 
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1. Introduction 



The instability of a self-gravitating incompressible fluid cylinder has been in- 
vestigated for first time by Chandrasekher and Fermi They have utilized 
the method of presenting the solinoidal vectors in terms of poloidal and toriodal 
quantities which is valid only for axisymmetric perturbations. Chandrasekher P] 
derived the stability criterion of such model for axisymmetric and non axisym- 
metric modes by using the normal mode analysis technique. Several extensions 
for such studies have been carried out upon considering the effect of different 
factors on the self-gravitating force 0, |^ . 

The purpose of the present work is to discuss the stability of a self-gravitating 
compressible fluid cylinder by utilizing the normal mode analysis. However, such 
technique is different from that used previously by Chandrasekher p[ and Radwan 
0], because the velocity field in a compressible fluid is no longer solinoidal, i.e. 



2. Basic equations 



We consider a selfgravitating fluid cylinder of radius embedded into a self 
gravitating tenuous medium of negligible motion. The fluid is assumed to be 
non-viscous and compressible. The model is acted upon by the self-gravitating 
and pressure gradient forces. We shall use the cylindrical coordinates (r, ip, z) 
with the z-axis coinciding with the axis of the cylinder. 

The basic equations appropriate for the present problem are the equations of 
motion, continuity equation, equation of state and the equations satisfying the 
self-gravitating potentials interior and exterior the fluid cylinder. These equations 
could be written in the form 



\jP + pyV, (1) 




d \ 

— + u-S7\p = -p{s7-u), (2) 



(3) 
(4) 




v' V' = 0, 



(5) 



where p, m, V and p are the fluid mass density, velocity vector, self-gravitating 
potential and kinetic pressure; is the self-gravitating potential of the tenuous 
medium exterior the fluid cylinder, 7 is the ratio of the specific heats and G is the 
gravitational constant. Here the superscript e denotes the exterior of the fluid 
cylinder. In the equilibrium state we have 



M =0, p = p , p = p , V = V V^ = V\ 



where index characterizes the equilibrium quantities, later on the index 1 is 
pertaining the perturbed quantities, equations (4) and (5) in this state read 



V; = -47rGp, (6) 



v'^/ = 0. (7) 

By integrating equations (6) and (7) with respect to r and determining the con- 
stants of integration. The latter may be determined upon applying the boundary 
condition that the self-gravitating potential and its derivative must be continuous 
across the boundary surface at r = R^^ . Consequently, the non-singular solutions 
of equations (6) and (7) are given by 



K = --^Gp r\ (8) 





r 



V; = -2TxGp RHn— + C (9) 



where is an arbitrary constant. In the present initial state, equation (1) yields 

VP =P VV (10) 

00" 

By integrating equation (10) and determining the integration constant (note that 
p = at r = i? ). Finally, the distribution of the pressure in the unperturbed 

state is given by 

p =nGp^iR^ -r^). (11) 

" 

It is worthwhile to mention here that p is not constant in contrast to other 



studies for different models acted upon capillary or/and electromagnetic forces, 
P], where it is found that p is constant. 



3. Perturbation Analysis 



For small departure from the equilibrium state, the fluid physical quantity Q 
could be expressed as 

g(r, z,t)^Yl e"^n(^, V:^:t), n = 0, 1 (12) 



where Q stands for p, u, V, and p. e is the amplitude of the perturbation at all 
times 

€ = € e"*, (13) 



where e (= e at t=0) is the initial amplitude and a is the growth rate. If a is 

^ 

imaginary say a — icu, {i — V—l), then — is the oscillation frequency. Based 

on the expansion (12) and from the view point of the linearized theory, the de- 
formation along the fluid cylinder interface due to perturbation could be written 
in the form 

r^R+R, \R \ « R , 

1 ' 1 ' 

with 

Here k and m are, respectively, the longitudinal and azimuthal wavenumbers and 
R is the elevation of the surface wave measured from the initial position at 
r = R. 



By the use of the expressions (12)-(14) for equations (l)-(5), the relevant 
perturbation equations are 

ap u = - ^p + p V K ' (15) 



<yp = -P (V ), (16) 

1 1 

IP 

p ^^p , (17) 
1 p 1 



V' K = 0' (IS) 
V^'V'^'^O. (19) 

By an appeal to ((/?, 2;)-dependence (cf. equation (14)) and based on the lin- 
ear perturbation technique concerning stability theory, every perturbed quantity 
Qi{r,ip, z,t) could be expressed as ei'^^+^C'^^+'^'p)) times an amplitude function of 
r. Consequently, the self-gravitating equations (18) and (19) are solved. Upon 
applying appropriate boundary conditions across the cylindrical fluid interface at 



A 



r = Ro, the constants of integration are determined and the non-singular solution 
of equations (18) and (19) are given by 



V = AnGp R K (x)I (kr)e 



(20) 



and 



V " = inGp R I ix)K (kr)e'^'''+"''f^ 



(21) 



Here / and K are the modified Bessel functions of the first and second kind of 

m m 

order m while x{— kR^) is the dimensionless longitudinal wavenumber. Combin- 
ing equations (15)-(17) and solving the resulting differential equations by utilizing 
similar steps as those which are used in solving equations (18) and (19). The con- 
stants of integration could be identified upon applying the kinematic boundary 
condition that the normal component of the velocity vector must be compatible 
with the velocity of the perturbed fluid interface. The non-singular solutions are 
given by 



1 



u — 

1 ap 



V 



AnGp'R Kjx)IJkr)-CIJrdk' + -) 







^i(,kz+m^) I ^ ^22) 



and 



where C is defined by 



Pi 



(23) 



P 



47rGp xK {x)r (x) - a' 



a 



a 



(24) 



and a(— 



IP 



-S-) is the speed of sound in the fiuid. 



Finally, we have to apply the boundary condition that the normal component 
of the stress tensor must be continuous across the fiuid cylindrical interface at 
r — R . This condition leads to 



dp 

p, + R ^ = 0, 
or 



at 



r = R 



(25) 



By substituting from equations (11), (23) and (24) into condition (25), following 
relation is obtained 



a 



^^Gp IJy) 



^Kjx)I'Jx)IJy) 1 



(26) 



where 



\ 



[a/R 



(27) 



is the compressible dimensionless longitudinal wavenumber. Here we get y = x 
as a — > oo. 

Equation (26) is the desired self-gravitating dispersion relation of an inviscid 
compressible fluid cylinder. It involves the entity {AnGp )~^/2 as a unit of time, 



the two kinds of the modified Bessel functions of different arguments, the sound 
speed in the fluid and the classical and compressible wavenumbers x and y. 

In the limiting case as a — oo while m > 0, the relation (26) reduces to 



AnGp- 



xl' (x) 



I (x)K (x) - 



{2i 



This relation coincides with the dispersion relation deduced by Chandrasekhar 
for solinoidal velocity vector. 

If we suppose that a oo and m = 0, the relation (26) reduces to 



AnGp 



xl^ {x) 
I ix) 



I (x)K (x) - - 



(29) 



The relation (29) has been derived for first time by Chandrasekhar and Fermi 
]![]. Indeed they have used a technique which is totally different from that used 
here. Such a technique is based on presenting the solinoidal vectors in terms of 
poloidal and toroidal quantities. 

As the compressibility factor influence is very small, the effect of the self- 
gravitating force on the instability of the fluid cylinder may be determined upon 
discussing the relation (27) and (28). 

By the aid of the numerical data of the modified Bessel functions [Q, it is 
found for all m 7^ that 

U^)Kjx)<i (30) 



Hence 



cr^ < for all m 7^ 0. 



(31) 



Therefore, the fluid cylinder is gravitationally stable for all purely non-axisymmetric 
perturbations. In the axisymmetric perturbation mode m = it is found that 



< 0, 



(32) 



in the domain < x < 1.0667 while 



a 



(33) 



at the critical point x — 1.0667 and 



a' > 0, (34) 

in the wide domain 1.0667 < x < oo. This means that the fluid cyhnder is 
self-gravitational unstable as long as the 

A > (27r/1.0667), (35) 

while it is stable iff the perturbed wavelength A satisfy 

A < (27r/1.0667), 

where the equality corresponds to the marginal stability state. 

Now, as we have seen in the foregoing discussion the self-gravitating incom- 
pressible fluid cylinder is unstable only in the axisymmetric mode m — 0. There- 
fore, in discussing the general relation (26) for determining the influence of the 
compressibility, we focus our stability discussions for in the axisymmetric pertur- 
bation mode m — 0. For m — 0, the relation (26) reduces to 



_yi[{y) 



xk^{x)r{x)I^{y) ^ 1 

yi'iy) 2 



(36) 



where the dimensionless compressible wavenumber y is still given by (27). The 
relation (36) has been calculated in the computer and the effect of different factors 
are identifled. Sec figure (1) and table (1). The analytical results are verified and 
it is found that the compressibility destabilizing the model for different values of 
a in particular as the sound speed a is very large. It is worthwhile to mention here 
also that for all values of s in the range 0.1 < s < 5.0 (with s — a/R^) we found 
that the model is unstable in the domain < x < 1.0233928. While it is stable 
in the neighboring domains 1.0233928 < ,x < oo where the equality corresponds 
to the marginal stable state. However, with increasing s values the areas under 
the unstable curves are increasing. Corresponding to s = 0.1,0.5,1.0,2.0 and 
5.0, the maximum modes of instability, respectively, are Nmax = 



0.066997,0.246781,0.309376,0.331429 and 0.338533 at x = 0.46,0.56,0.61,0.61 
and 0.66. This means that the compressibility has a tendency of destabilizing 
the model. On contrary as s = we have the selfgravitating unstable domain 
is < X < 1.0678. This shows that all the unstable domain is the presence 
of compressibility are less than that in the absence of the compressibility. This 
means that the compressibility has a tendency of stabilizing the model. 
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Table 1 



Stable and unstable domain for selfgravitating compressible fluid cylinder with 
S — — and UJ = . 

^0 



n 



Figure (1) 



Stable and unstable domains for selfgravitating compressible cylinder with 

—, , and cu* = , 
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